Impurity band in clean superconducting weak links 
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Weak impurity scattering produces a narrow band with a finite density of states near the phase 
difference = tt in the mid-gap energy spectrum of a macroscopic superconducting weak link. The 
equivalent distribution of transmission coefficients of various conducting quantum channels is found. 
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Effect of impurities on transport properties of super- 
conducting nanostructures is one of the key issues in 
the physics of mesoscopic superconductors (see 0] for 
a review). Recent advances in nanotechnology revived 
interest in quantum point contacts (see |^ and refer- 
ences therein) and in other devices that employ quantum 
conductors connected to superconducting electrodes d). 
SNS junctions and weak links consisting of two supercon- 
ductors connected by a small orifice in a thin insulating 
layer (point contacts) are the simplest devices of interest. 
Though impurity effects in SNS junctions are well stud- 
ied (the work began with Ref. and still goes on, see 
for example ||]), the role of impurity scattering in super- 
conducting point contacts remains not fully investigated. 

As an example, consider a ballistic point contact be- 
tween two superconductors assuming that the thickness 
d of the insulating layer and the size of the orifice a are 
shorter than the coherence length ^ and the impurity 
mean free path £. As is well known Q there exist 
mid-gap states with the spectrum 



±eo, en = |A|cos(0/2) 
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Here (j) = X2 ~ Xii with xi,2 being the order parameter 
phase on the left (right) from the orifice; the upper sign 
refers to particles moving to the right from region 1 into 
region 2 and vice versa. Impurities do not appear in Eq. 
(|l|) because the characteristic dimension is determined by 
the size of the constriction d rather than by the param- 
eter d/i being assumed infinitely small. It is natural to 
expect however that the impurity scattering would mod- 
ify this spectrum at an energy scale determined by the 
small parameter d/i, especially for low energies where the 
two branches of the spectrum for right- and left-moving 
particles cross. The same can be expected for a long bal- 
listic SNS junction which has a mid-gap spectrum |Q, |^ 
consisting of many branches for right- and left-moving 
particles that cross at (f) — ir. 

In this Letter we consider both point contacts and long 
SNS junctions and show that a weak impurity scattering 
transforms their mid-gap spectra in such a way that a 
narrow band having a finite density of states and a width 
~ |A|-y/ d/£ appears near each crossing point at the phase 
difference (j) — n. This impurity band is expected to 
have a dramatic effect on dynamic properties of weak 



links. In particular, it enhances the inelastic electron- 
phonon relaxation rate at low temperatures which would 
be otherwise nearly zero due to the energy conservation. 

Impurity band in a point contact.-We consider first a 
ballistic point contact such that a ^ d <^ £^ <^ £. One 
has Ai.2 = lAje*-^"^'^ to the left (right) from the orifice, 
respectively (see Fig. |l| ). The quasiclassical Green func- 
tions (retarded or advanced) satisfy the normalization 
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ff^ = 1 and obey the Eilcnberger equations 
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where s is the distance along the particle trajectory. We 
assume zero magnetic field. The right-hand sides of these 
equations describe the scattering by impurities. We use 
(...) to denote an average over the Fermi surface. The 
standard technique of averaging over impurities is appli- 
cable because the number of impurities within the vol- 
ume of the orifice is large. Indeed, the mean free time is 
r^i ^ A^F'T^impl^P where \u\ ~ Upp^ is the Fourier trans- 
form of the impurity potential U , and Np is the normal 
single-spin density of states at the Fermi level. The num- 
ber of impurities nimprfa^ ~ {d/£){ppa)'^{Ep /U)"^ can be 
very large even for U ~ Ep because of a macroscopic 
number of quantum channels in the orifice {ppaf' ^ 1. 

For low energies |e| ^ |A| and close to = tt the 
Green functions are localized near the orifice at distances 
s ~ ^. We put = TT -f (5 where (5 ^ 1 so that Ai = 
-z|A|e-*''/2, A2 = i|A|e*''/2. Following § we write 



(5) 



The upper sign is for right-moving particles, px > 0, the 
lower sign is for left-moving particles, Px < 0. The x 
axis is perpendicular to the insulating layer. We assume 
^ 1 — 77^ so that g± = i(±. Equations (||, ||) become 

vp^ + 2C±|A| cos((5/2)sign(s) = 0, 

vp^ - (2e T 2|A| sm{d/2) + ia±) C± = 0, 
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FIG. 1: The point contact between two superconductors (Si 
and S2) separated by an insulating layer (I). The x axis is 
directed from the region where the order parameter phase is 
Xi into the region where the phase is xa- 



where 



-i((5>±(/>/2±(/t)/2) 
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The boundary conditions for s — > cxd follow from the 
expressions for the Green functions in the bulk. For small 
(5, the functions C± ^-nd ?7± — s- 1. The solution is 
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C±expi-K), K^2\A\\s\/vF, 

C±vp^ [ e-^ {2\A\E± + ia±) ds, 
Jo 
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E± + ivp^ / a± ds 
Jo 



Here E± = e/\A\ T sm{S/2). 

The averages (g), etc., are proportional to the solid an- 
gle at which the orifice is visible from the position point, 
they decrease quickly at distances s ~ a from the orifice. 
Indeed, g (p, r) = =f/ (p, r) — i( (0) for trajectories that 
go through the orifice and g (p, r) = / (p, r) = oth- 
erwise because, for non-through trajectories, the Green 
functions are small as compared to ^ 1. For s ^ a one 
has K = 0. Moreover, using Eq. (0) we find 



a±ds — £ 



ds(g(0))^=75T(0) 



where = {2tt) ^ J^^^^^g{p,r) dilp, and 7 - d/e is 
a geometric factor depending on the shape of the contact 
and on the position of the trajectory with respect to the 
orifice. We neglect the latter dependence and consider a 
constant 7 for simplicity. Finally, 
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In the limit 7 — s- 0, the Green functions have poles at 
E± = which is Eq. (]l|). The spectrum is shown in Fig. 
H by lines (1) for > and (2) for < 0. However, as 
E+E- approaches 7, the apparent pole-like behavior of 
Eq. transforms into a more complicated dependence. 
To calculate the Green functions, we solve Eq. for 
5±(0). If < E+E^ < 47, i.e.. 



sin2(V2) < eViAp < 47 sm'^{S/2) 



(8) 




FIG. 2: Lines 1, 2, 3, and 4 encompass the energy band 
with a finite density of states (shaded region). The maximum 
half-width of the band is 2^/7. For a ballistic contact r — > 00, 
the spectrum is given by Eq. (^ and follows lines 1 and 2. 
The entire pattern is 27r-periodic in 0. 



we find g^! (s) = g^ (0) e ^ where 



5^(0) 
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and g^ — —[g^]*. The radical is defined as an analytical 
function with the cuts along the borders of the region 
determined by Eq. (||) (shaded region in Fig. ||). The sign 
is chosen for e within the upper part of the region. The 
normalized density of states Ve{s) = [g^ (s) — g^ (s)] /2 
is nonzero within the energy interval of Eq. (^. The 
maximum half- width of the energy band 2^/j is reached 
for (5 = 0. Inside the constriction 



(10) 



For a particle with a given sign of p^ the density of states 
is nonzero near both E- = and E+ — 0: the scattering 
mixes states with positive and negative px- However, far 
from the crossing point, |e| ^ \A\y/j, the ratio i'±/i'zf: = 
E^/E± is large near the corresponding ballistic spectrum 
E± = 0: its magnitude is of the order e^/|Ap7. 

Beyond the range of Eq. (H) g^ and g^ coincide and 
the density of states vanishes. For E+E^ > 4j, 
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For 7 <C \Ej^E- \ one recovers the poles g^'^^ (0) — i/E± 
with the density of states v = 7r|A|(5(e =F eo)- Obvi- 
ously, Eq. ( p^ ) conserves the total number of states since 
J v±de — 7r|A| for any (p. 

Long SNS bridge - Consider now a normal bridge of a 
length d ^ that connects two massive superconductors. 
Its width is much shorter than ^. We assume that it has 
specular walls. Irregularities on the walls can also be 
modelled by random impurities. Superconductors and 
the normal metal have the same Fermi velocities vp- We 
take the x axis along the bridge (see Fig. ||). 



FIG. 3: Normal bridge between two bulk superconductors. 



Let s — zLsq — id/I cos6'| be the outlets of the normal 
bridge into the bulk superconductors. For s < —sq the 
solutions in the bulk are ^ 

9±=9oo+ .g±fee'=(^+^«), /± = /oo + 5±fc/3fee'=(^+^«) (11) 

where = A/ty/\A \^ - g^o = e/VI^P " ^hile 

k = 2?;^^y^|A|^^~^ and /3fc = A/ {e + ivFk/2). The 
impurity scattering in the bulk can be neglected. The 
signs correspond to ±px > 0; the order parameter phase 
should be xi for Px > and X2 for Px < 0. For s > sq 



9± = 5o 



9±ke 



, f± ^ + 9±kf3-ke-'^'^'«l 

(12) 



Here the phase is X2 for Px > and xi for < 0. 

We first outline the known solution ^ for states 
with not very low energies vp/d <^ e such that d ^ 
Here ^at ~ vpjT. In the region inside the bridge where 
A = Eqs. (||-^ yield = g]^ ^ const and 



/ = C exp 



2e 
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where Vx — vpcosO and x — scosO. We assumed that 
(./) = (/O ~ Indeed, / and oscillate rapidly as 
functions of 9 since ed/vp 3> 1 and vanish after averaging 
over the angles. The continuity at the borders between 
the bridge and the bulk gives g±k = g±N - 

tan a ± igoo 



and 



9±N 



goo tan a ±i 



(13) 



where a = [d/vx) {'2e + i {gn) /t) - 4>/2. For |e| < |A|, 
the function g±M has poles when a — ± arccos (e/ | A|) + 
Trn or I, § 
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where uJx — \vx\/2d. The state with n = —1 has e = for 
(j) — n. The spectrum for t = oo is shown in Fig. ^ by 
lines 1 and 2. For e ^ vp/d, the poles are closely packed. 
As a result, the angular average collects contributions 
from many poles which givesQ] ^3^'"^'^ = ±1. 

For low energies e ^ |A|, the term goo tana <C 1 if d » 



^. Eq. (13) results in g^ = tana. The Green function 



FIG. 4: The energy bands (shaded regions) for a long SNS 
bridge. Lines 1 and 2 show the spectrum of a particle with 
Pa; > and Px < 0, respectively, in a ballistic bridge £ oo. 
The pattern is 27r-periodic in (j). 



is large when a is close to 7r/2 + irn, i.e., close to the 
points e — uonm where the branches for px > and for 
Pa: < cross. Cousidcr these regions in more detail. We 
put e = e' + LUomn and write the functions in the normal 
bridge as / = f ex.])(ztiTrmx / d) , = exp(=pz7rma;/d). 
We now express / and /t through rj and according to 
Eq. (^ and assume gpf = iQ. The solution of Eqs. (||, ^) 
is Q± = A± ~ const and 



(2e' + ia±) dx. 



Here a± is determined by Eq. (^ where /, are replaced 
with /, p. Note that now K — Q. Matching this with 
Eqs. (O, O) aX X ^ ±d we obtain 



9±N 



(14) 



e' + i(T±/2 =F t^x5/2 

where (j) = -k + 5. Using Eq. (||) we find a± — (g)-^ jr. 

The angle-resolved density of states is proportional to 
the averaged g^' — g^ ■ Calculating the average we obtain 



{gN)± = «wo 



/i d[L 
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where cjo = vpl2d. We put x — 5/2^, y — 
e'/ijJoV7i -^i = \/i{9n)± where 7 = d/i'. Performing 
integration in Eq. (|l^) for 5 ~ ^77 and e ^ ^0^^ we find 



(16) 



For small a; one has z+ = z_ = ± \J2 — ip /2 

which is similar Eq. (^. For y < ^/2, the functions 
(5^) = ^ {9^) are complex, thus the density of states 
is nonzero, {v) ~ l/^/j. For y > solutions are imag- 
inary, V = 0. To find the limit where a nonzero density 
of states appears, we look for z = ia and find conditions 
when there are no real solutions for a. We have 
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For small x we return to the previous result: the density 
of states vanishes for y > \/2. However, if e' = y = a 
real solution for a exists for any x: the density of states 
is zero along the axis e' = which agrees with Eq. (H). 
Next, we note that the density of states is singular at one 
of the borders of its region of existence. Thus, either a+ 
or a_ should go to infinity. Assume that a-|_ —^ oo. Then 
a_ — > according to the second equation (pT[). In turn, 



Eqs. ® and (|T|) we find 



= ^ In 

X x'^ 



y 



Infinite values for a-|_ are possible for y = x. Similarly, 
assuming a_ — s- oo, a_|_ — > we finally find the borders 
of the region with a nonzero density of states, y — or 



uo\5\/2. 



(18) 



Equation (|l|) gives the lower-|e'| limits shown in Fig. M 
by lines 1 and 2. The higher-energy limits (lines 3 and 
4 in Fig. ^ are set by \y\ = \/2 for small x. For large x 
and y we expect that one of a± is large while the other is 



small. We then again obtain Eq. (18) as the asymptotic 
expression. Therefore, the two borders approach each 
other for large x and y, as was also the case for a point 
contact. We thus expect that Eq. ([lo| ) provides a qual- 
itatively correct representation for the averaged density 
of states of a long SNS junction, as well. 

Discussion.-ThB results for point contacts and for long 
SNS bridges are qualitatively similar. To simplify the dis- 
cussion we consider the impurity band in a point contact 
where an explicit equation for the density of states is 
available. It can be compared with the known mid-gap 
energy spectrum0, |ll| for a contact that has a tunnel 
barrier with a transmission probability Z?, 



e = ±ec, (D = |A|a/1 ~Dsm\(jy/2) 



(19) 



In Fig. gthis spectrum follows lines (3) and (4) and has a 
gap with the half- width | A|\/i? at = tt where R—l — D 
is the reflection coefficient. The gap results from the mix- 
ing of right- and left-moving particles provided by the 
barrier. One can consider the impurity band Eq. (^ as 
a result of superposition of a large number of conduction 
channels with various transmission coefficients D result- 
ing from scattering of particles on different trajectories 
characterized by a given momentum. Superposition of 
states with different momenta due to scattering by im- 
purities is also known to produce an impurity band in 
d-wave superconductors near the gap nodes [l^ . To find 
the equivalent probability distribution P{D) we write the 
total density of states for an energy e in the form 

+ {i^)_=tt\A\J [S{e~eD)+S{e + eo)] P{D) dD. 

The relevant values of D are close to one. With help of 



1 a/4'7 - R 
P{D) = —^^^L^, i?< 47. (20) 
27r7 Vi? 



The distribution is truncated at i? > 47. For 7 — > 
one obtains P{D) — S{R). The square-root singularity 
at i? = in Eq. (20) resembles that of the universal 
distribution [ p^ for diffusive normal conductors. The 
singularity ensures a finite density of states at (p — tt 
which is J/ = — (e^/47|AP) according to Eq. 

(0). However, Eq. (f^) is not universal: in contrast to 
its diffusive counterpart it depends on the geometry of 
the contact and on the quasiparticle mean free path. 

The impurity band is expected to have a profound 
effect on dynamics of weak links. For example, it en- 
hances the electron-phonon relaxation at low tempera- 
tures which otherwise would almost vanish due to the 
energy conservation. Indeed, the electron-phonon relax- 
ation rate is proportional to the product of two electron 
densities of states and i^ei, and the phonon density of 
states iVph(w) where uj = e — ei. Since A^ph(w) cx the 
product vanishes if each electronic density of states is a 
delta function S{e— cq) with only one state eo for a given 
phase (p and the sign of momentum . It is the impurity 
scattering that broadens the density of states into a band 
thus making the electron-phonon relaxation possible. 
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